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obtain a surprisingly predictive scenario. The absolute neutrino mass scale, the Dirac and the two 
Majorana phases of the neutrino mixing matrix remain determined in terms of the set of already 
measured low energy observables, modulo a discrete ambiguity in the signs of two neutrino mixing 
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I , I. INTRODUCTION 

<D ■ 

In the last decade, experiments with solar, atmospheric, reactor and accelerator neutrinos have provided compelling 

evidences for neutrino oscillations [l|, which imply nonvanishing neutrino masses and mixings. Neutrino oscillation 

experiments have been measuring with increasing precision the values of the mixing angles and of the mass-squared 

differences. The value of the absolute neutrino mass scale is still unknown; however, existing limits imply that this 

scale is baffiingly small, much smaller than those of all the other elementary fermions. The most popular explanation 

r»o ■ for the neutrino mass suppression is undoubtedly provided by the see-saw mechanism [2| which requires the existence 
of very heavy right-handed (RH) Majorana neutrinos. Fermions with quantum numbers of RH neutrinos, that are 
singlets under the standard model (SM) gauge group, are found in the spinorial 16 representation of SO(10) [j| |J], 

^f-\ which therefore provides a quite natural Grand Unified Theory (GUT) framework to embed the see-saw. 

The see-saw RH neutrinos also play a key role in leptogenesis [5|, l|3| , which is a very appealing scenario to explain the 
origin of the Baryon Asymmetry of the Universe (BAU). In leptogenesis, the cosmological baryon asymmetry is seeded 
by an initial asymmetry in lepton number generated in the out-of-equilibrium decay of the RH neutrinos, that is then 
transferred in part to baryons by means of the B + L violating 'sphaleron' interactions, that are non-perturbative SM 
processes. In SO(10), the order of magnitude of the RH neutrino masses is fixed around the scale of the spontaneous 
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breaking of the B — L U(l) symmetry, and it is consistent with the values of the RH neutrino masses required for 
successful leptogenesis Mr ~ 10 11±2 GeV. Indeed, the double role of RH neutrinos in the see-saw and in leptogenesis 
underlines the importance of deriving information on their mass spectrum. 

Recently, an analysis of the relations between the left-handed (LH) neutrino observables (mass-squared differences 
and mixings) and the RH neutrino spectrum, constrained to be of a compact form (i.e. with masses all of the 
same order of magnitude) was carried out within the framework of an SO(10)-inspired model [7|, and a scenario 
for baryogenesis via leptogenesis was also constructed. The study in Ref. [7j was carried out under the simplifying 
assumption of a vanishing value of the lepton mixing angle #13, and in the leptogenesis analysis all lepton flavour 
effects [a-ROll as well as the effects from the heavier RH neutrinos [IJ, [l2j had been neglected. However, recent 
experimental results hint to a nonvanishing value of #13 [13l - ll5| and imply that the assumption #13 = should be 
dropped. The inclusion of flavour effects is also mandatory when leptogenesis occurs below T ~ 10 12 GeV, since the 
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one flavour 'approximation' is known to give unreliable results. Moreover, in the case of a compact RH spectrum, 
that is when all the RH neutrino masses fall within a factor of a few, to obtain a trustworthy result it is also necessary 
to include the asymmetry production and washouts from the two heavier RH neutrinos. 

In the present paper we consider a scenario similar to the one in Ref. [7]] improving on several points. We fix #13 
to the nonvanishing best fit value given in Ref. [15j . This in turn implies that the Dirac phase 8 of the Pontecorvo, 
Maki, Nagakawa and Sakata (PMNS) mixing matrix [111 LLZ| enters all the equations, and in particular contributes to 
the leptogenesis CP asymmetries. Most importantly, we clarify how the conditions ensuring a compact RH neutrino 
spectrum have consistent solutions only for <5 ^ 0, and how the corresponding solutions yield a surprisingly predictive 
scenario in which all the yet unknown low energy parameters, namely the LH neutrino mass scale mi and the three 
PMNS CP violating phases 5, a and (3, remain determined in terms of already measured quantities, modulo a few 
signs ambiguities. In the high energy sector, the RH neutrino spectrum is also predicted. The crucial test of the 
scenario is then the computation of the baryon asymmetry yield of leptogenesis. We include lepton flavour effects 
|8l4l0l| in our analysis and argue that they are crucial to evaluate correctly the baryon asymmetry. Most importantly, 
the high level of predictability of our framework allows to predict both the size and the sign of the BAU. By requiring 
agreement with observations, we are then able to solve almost completely the residual signs ambiguities. 

The paper is organized as follows. In Section [II] we describe the SO(10) framework and spell out the quark-lepton 
symmetry assumption. In Section [IIII we discuss the constraining conditions that ensure a compact spectrum for the 
RH neutrinos. In spite of quark-lepton symmetry a compact form is achieved, and with a sufficiently large scale not 
to conflict with the Davidson-Ibarra bound that, within the 5*0(10) see-saw, often vetoes successful baryogenesis via 
leptogenisis. In Section IIVI we confront our scenario with the set of measured low energy neutrino observables, and 
we work out predictions for the absolute scale of neutrino masses mi, for the PMNS CP violating phases S, a, (3 and 
for the RH neutrino mass matrix. In Section [V] we calculate the various CP asymmetries in RH neutrino decays, we 
briefly discuss the procedure followed to estimate the baryon asymmetry yield of leptogenesis and stress how a proper 
treatment of flavour effects is crucial for obtaining reliable estimates for the different cases. Finally in Section IVll we 
discuss our results and draw the conclusions. 



II. SO(10) GUT AND QUARK-LEPTON SYMMETRY 

We work in a non-supersymmetric grand unified SO(10) model. We assume three fermion families whose left- 
handed (LH) states are assigned to a 16 spinorial representation of SO(10), which thus contains all the SM fermion 
and antifcrmion states of the same chirality. In addition, the 16 includes one SU{2) singlet neutrino for each family. 
All elementary fermions of opposite chirality (RH) are assigned to the conjugate representation 16. Fermion masses 
are generated by Yukawa terms of the form: 

Y i:i ■ 16, H^ 16j + h.c. (1) 

where Y is a 3 x 3 matrix of Yukawa couplings with indexes in family space, and H denotes a multiplet of scalar 
(Higgs) bosons. The tensor product of the fermion representations in eq. ([1]) gives 

16 x 16 = 10 s + 126, + 120 Q (2) 

where the subscripts s, a refer to the symmetric and antisymmetric nature of the representation in the family indexes. 
Thus, to make the Yukawa term in eq. (HJ an SO(10) singlet, H must be assigned to a 10, to a 126, or to a 120. 
Clearly, for 10 and 126, that match the two symmetric fragments of the tensor product, the Yukawa matrix Y is 
symmetric, while for the 120 it is antisymmetric. 

In the present work we consider Yukawa terms involving only the 10 and 126 that are already needed for the gauge 
symmetry breaking pattern (the 126 is also needed to generate Majorana masses for the see-saw RH neutrinos) and 
we exclude Yukawa couplings with the 120 which would imply a departure from minimality. 

We adopt the following SO(10) breaking pattern; we also indicate the set of Higgs SO(10) representations needed 
for each step 

5O(10) ^> SU(2) L <8> SU(2) r <g> 5*7(4) 
i^> SU(2) L x U(1) Y x SU(3) C 
12 °^ SU(3) C x U(1) Q (3) 



The first step in this chain is a breaking to a maximal subgroup of 5O(10), the intermediate Pati-Salam group 



SU(2) L ® SU(2) R ® SU{4)^, Let us list explicitly useful branching rules for SO(10) -> SU(2) L ® SU(2) R ® SU{4): 

210 D (1,1,1)8... 

126 = (1,1, 6) 8 (3, 1,10) ©(1,3, 10) ©(2, 2, 15) 

16 = (2, 1,4) ©(1,2, 4) 

16® 16 D (1,3, TO)©... (4) 

The SU(2) L ® SU{2) R <g> SU(4) singlet in the 210 is responsible for the first breaking at the GUT scale. 

The 126 then breaks SU(2) L ®SU(2) R ® 517(4) -> SU{2) L ®U{1) Y ® SU(3) C at an intermediate scale A R , that we 
assume to be around 10 11 GeV. The relevant component that triggers the breaking is (1, 3, 10) since the 10 of 5[/(4) 
contains an 5C/(3) singlet. The other components, (3, 1, 10) and (2, 2, 15) would break SU(2)l, while (1, 1, 6) would 
break colour. The RH neutrino N, together with all the other SM SU(2)^ singlet fields, is contained in (1, 2, 4), and 
thus the bilinear TV • N belongs to the fragment displayed in the last line, which is the only one suited to build up 
a gauge invariant term when coupled to the intermediate gauge symmetry breaking component (1, 3, 10) of the 126 
Higgs multiplet. We can also decompose 50(10) according to 50(10) D 517(5) ® 17(1) D 517(5). With respect to 
5C/(5): 

126 = 1 + 5 + 10 + 15 + 45 + 50. (5) 

Of these representations only the 1, 5, and 45 have neutral colour singlet Higgs components that can have nonzero 
vacuum expectation values if U(l) em x SU(3)c has to remain unbroken. N is a singlet with respect to 5/7(5) and so 
it is N ■ N; therefore, it is the SU(5) singlet in 126 that couples to N ■ N and gives an invariant mass of the order of 
the intermediate scale to N. In terms of the intermediate representations, it is (1, 3, 10) in last line of Eq. (|4|, that 
must contain an SU(5) singlet. 

As regards the fermion masses, if they originate only from vacuum expectation values (vevs) of scalars in the 10, 
the following relations hold: 

mo = m u , m e = rrid, (6) 

where mo is the neutrino Dirac mass matrix, and m u , rrid and m e are respectively the mass matrices for the up and 
down quarks and charged leptons. The two relations in eq. ([6]) are sometimes referred to as quark-lepton symmetry; 
they imply for each generation the GUT scale prediction rriejmdi = 1 (with i = 1,2,3 a generation index) which 
however, after including renormalization group corrections, agrees with observations only for the third generation (b-r 
unification) but is badly violated for the first and second generations. If instead quark and lepton masses originate 
from one or more 126 (that however should be different from the 126 that breaks the gauge symmetry) the following 
relations hold: 

m D = -3ra„, m e = -3m d . (7) 

The factor of —3 is a colour factor between leptons and quarks; it is reminiscent of the Georgi and Jarlskog mecha- 
nism [19], which in SU(5), when it appears in a family dependent way, allows to circumvent the prediction of unification 
for the first two families yielding m^/rris ^ m e /md ^ 1 while preserving m T /mb — 1. In SU(5) the discrepancy with 
the observed values of the down-quarks and charged lepton masses can in fact be weakened by assuming that the 
Yukawa coupling of the second generation to itself involves a 45 of SU(5) scalars, instead of the usual 5, yielding mass 
ratios 

|"fy/ m «l = l"W TO el = 3 ( 8 ) 

which are in better agreement with the measured values. Now, under SO(10) — > SU(5) x U(l) the 126 contains 
precisely a 45 which, as Harvey, Ramond and Reiss showed [20, [21j , allows to implement the same mechanism also 
in SO(10). 

In our SO(10) model we assume that the SM SU(2)l Higgs doublet is a combination of representations in the 10 
and 126 of SO(10). While this allows to account for non-unification for the down-quark and charged lepton masses 
of the two lightest families, it still predicts an approximate quark-lepton symmetry, that is m e ~ mj are in any case 
connected by coefficients of order 1. As regards the neutrino Dirac mass matrix, for definiteness we will stick to the 
simpler relation in eq. (j6|), which can coexist with eq. (0 for m MiS , m e ,d if the u-v sector masses are dominated by 
the 10 vevs. We also assume that, in the diagonal basis for the down-quarks and charged leptons mass matrices, 
the unitary rotation Vl that diagonalizes the symmetric matrix nip coincides with the Cabibbo-Kobayashi-Maskawa 
(CKM) rotation that diagonalizes m u . Namely we assume as a working hypothesis: 

m D = m u and V L = V C km ■ (9) 



We stress at this point that our results do not depend in any crucial way on the precise form of the quark-lepton 
relations, and the ansatz eq. ([§]) is adopted here only for the sake of simplicity. However, the possibility of constructing 
a predictive framework does depend on the fact that in SO(10) a precise relation between mo and m u and Vl and 
Vckm exists, which naturally follows from fermions unification within a single irreducible representation of the group. 
In fact, once the details of the symmetry breaking pattern and of the fermion couplings to the 10 and 126 are given, 
a quark-lepton mass relation remains in any case fixed, and in particular a highly hierarchical spectrum for the 
eigenvalues of mn is a straightforward consequence of the SO(10) GUT framework (see also Ref. |22|]). As regards 
the full 6x6 mass matrix of the neutral sector, recalling that symmetric Yukawa matrices imply that mjj = mo, it 
can be written as 

M=( ° £»V (10) 

y m D M R J y ' 

where mo and Mr receive respectively contributions from the following vevs: 

m D ~ (10 + 126) Aew , M R ~ (126) Ar , (11) 

where the 126 contributing to Mr has a vev O(Ar) 3> Aew along the 517(5) singlet component. 

III. COMPACT RH NEUTRINO SPECTRUM 

The hierarchy Aew/Ar <C 1 between the two types of vevs in eq. (|11[) enforces the see-saw mechanism, and, after 
diagonalizing the matrix (I10[) . one obtains the light neutrino mass matrix m„ from the seesaw formula that, with 
mjj = mo, reads: 

m v ~ ~mo M^ 1 mo ■ (12) 

Inverting the seesaw formula (fT2"|) gives 

Mr ~ —mD m" 1 -hid, (13) 

which shows that one can obtain information on Mr by using the available experimental data on m„, and assuming 
quark-lepton symmetry for mo- 

Quark-lepton symmetry however, rendersproblematic the implementation of the mechanism of baryogenesis via 
leptogenesis within the SO(10) see-saw [23, [24J. This is due to the two factors of mo in eq. (|13l) that in general yield 
a very hierarchical spectrum for the RH neutrinos. In fact, by fixing the intermediate scale Ar around 10 11 GeV, the 
lightest RH state iVi , which is generally the main one responsible for generating a lepton asymmetry, acquires a mass 
Mri "C 10 9 GeV, that is well below the Davidson-Ibarra (DI) limit [25( which gives the benchmark to guarantee a 
sufficient production of lepton asymmetry from RH neutrino decays. There are basically two ways out to this problem. 
The first one relies on the fact that under certain conditions leptogenesis can also proceed via the decays of the two 



heavier RH neutrinos 11], whose masses remain well above the DI bound. Refs. [26143 1| present specific realizations of 



this possibility. The second way out relies on the possibility of enhancing resonantly the CP asymmetries [32j , which 
allows to evade completely the DI bound, but requires that at least one pair of RH neutrinos is highly degenerate in 
mass. In this paper we explore a third possibility, namely that in spite of the quark-lepton symmetry, the RH neutrino 
spectrum could still turn out to be of a compact form [7] that is, characterized by at most mildly hierarchical mass 
eigenvalues, all with values within the range 10 11 GeV which is the optimal one for leptogenesis. Clearly, such a 
possibility would avoid from the start the problem of a too light N\ . Let us see in detail how this possibility can be 
implemented. 

A. Conditions for a compact Ni spectrum 

A generic Dirac neutrino mass mo can be diagonalized by means of a biunitary transformation with two unitary 
matrices Vl and Vr, that is mo — V^jti^Vr. However, due to the assumed specific symmetry breaking pattern, in 
our model mo is symmetric, and in this case there exists (Takagi factorization 133]) a single unitary matrix Vl such 
that 

m D = v£m d j"Vl, (14) 



where m 



diag 



diag (mDi,TriD2,TnD3) is diagonal with real and non-negative eigenvalues. It follows that the RH 



neutrino mass matrix can be written in the form 

Mr = 

where we have introduced the symmetric matrix 



■Vlm% aa Am d ™ 9 Vl 



A = V??n- 1 v} 



(15) 



(16) 



On a naturalness ground, from the current knowledge about the light neutrino mass matrix m„, in the basis where 
charged lepton mass matrix is diagonal one would expect that the elements of A are at most mildly hierarchical. 



Then, if 



diaq 
l D 



is hierarchical as implied by quark-lepton symmetry, we would generally obtain a hierarchical RH 



neutrino spectrum. Therefore only a quite specific structure of the A matrix in eq. (|16p can enforce the conditions 
that ensure that the RH neutrino spectrum is compact. To illustrate this issue, let us recall that we are working under 
the assumption of quark-lepton symmetry eq. ([9]) which implies in particular that in the basis where the mass matrix 
of the charged leptons and of the down- type quarks are diagonal Vl = Vqkm ■ Although the results in Sections IIVI 
and fVl are obtained with the assumption Vl = Vckai, to write down reasonable analytical expressions for the RH 
neutrino mass spectrum in this section, we will set in first approximation Vl = ^3x3 (where the ^3x3 is the 3x3 
identity matrix). Eq. (fT5|) then gives 



(A n m 2 D1 A 12 m D1 m D2 Ai 3 m D1 m D3 
A 12 m D1 m D2 A 22 m 2 D2 A 23 m D2 m D3 
A 13 m D1 m D3 A 23 m D2 m D3 A 33 m 2 m 



(17) 



with A s=a m v x (recall that because A is symmetric A^ — Aji). Quark-lepton symmetry implies to_d3 3> xnoi 3> rnni, 
which suggests that a generically compact RH spectrum would result if 
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m D2 
m D3 



(18) 



since, if this were the case, all the hierarchically large entries in Mr would be sufficiently suppressed (notice that 
the Mr(1,3) and Mr(2,2) elements of eq. (fTT|) are non-hierarchical because from the light neutrino mass matrix we 



expect A 22 ~ Ai 3 and from the quark-lepton symmetry we expect mDimos 
assume the more restrictive condition: 



i D2 ). However, in this paper we will 
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« 1, (19) 



since they are needed to justify the simplifying approximation in eq. (|20[) below. The interest in exploring a scenario 
in which the two conditions in eq. (Tl9t are realized stems from the fact that it is quite likely that an A-matrix of this 
form would render leptogenesis a viable mechanism to explain the BAU within the SO(10) seesaw framework. In this 
paper, we will not speculate on the possible origin of the two relations in eq. (fT5|) . nor we will attempt to reproduce 
them by starting from a suitable fundamental Lagrangian, and thus the possibility that such a pattern could arise 
basically relies only on the fact that similar hierarchies do exist among quantities related to the Yukawa coupling 
sector. In fact, we believe that building up a theoretical justification for A 33 /An, A 23 /An <C 1 could be equally 
difficult than explaining the mass hierarchies of the charged fermions, which is a long standing unsolved problem in 
particle physics. We will, however, prove that as long as #13 is nonvanishing, imposing such relations is a technically 
consistent procedure, in the sense that they can be always fulfilled, regardless of the specific types of quark-lepton 
symmetry relations assumed. If A 23 and A 33 are negligible with respect to all the other entries, we can set in first 
approximation 

A 23 = A 33 =0. (20) 

As we will see, from these two conditions it follows that, besides obtaining a compact RH spectrum, two eigenvalues 
in Mr will actually be close to degenerate. In general, the degeneracy of pairs of Mr eigenvalues represents an 
interesting situation for leptogenesis, since it can allow for resonant enhancement of the CP asymmetries. Although 
we will find that, eventually, conditions eq. (|2T)|) are not sufficient to bring the dynamics of leptogensis fully within 
the resonant regime, it is still worth studying which class of general conditions for A could yield a pair of eigenvalues 
very close in mass 1 . 



1 The fact that our results for the Cosmic baryon asymmetry do not benefit from resonant enhancements of the CP asymmetries justifies 
the claim that the compact RH spectrum scenario represents a third possibility for realizing leptogenesis within the SO(10) GUT. 



The eigenvalues A of Mr in eq. (fl7|) are given by the solutions to the characteristic cubic equation 

X 3 + b\ 2 +c\ + d = 0, (21) 

with 

b = A n m 2 D1 + A 22 m 2 D2 + A 33 m 2 D3 , 

c = (A n A 2 2 - A\ 2 ) m 2 D1 m 2 D2 + (A n A 33 - A\ 3 ) m 2 D1 m 2 D3 + (A 22 A 33 - A\ 3 ) m 2 D2 m 2 D3 , 

d = (2A 12 A 13 A 23 + A n A 22 A 33 - A 2 13 A 22 - A 2 12 A 33 - A 2 23 A n ) m 2 m m 2 D2 m 2 D3 . (22) 

The necessary condition for two eigenvalues being equal is that the discriminant of eq. (|2ip vanishes. We can write 
down the discriminant as follows 

A = b 2 c 2 - 4c 3 - 4b 3 d + 18bcd - 27d 2 

« (^23 - A 2 2A 33 ) 2 m% 2 m 4 m [A 2 22 m% 2 + A 2 33 m% 3 + 2 (2A 2 23 - A 22 A 33 ) m 2 D2 m 2 D3 ] , (23) 

where in the second line we have expanded up to first order in "!^sx ^ ShL We have A = if 

r f mD3 mt 

We will consider only the first possibility, that involves solely elements of the matrix A. We then see that if A 23 , A 33 ps 
or alternatively A 23 , A 22 ps 0, quasi degeneracy of two RH neutrino masses results. Notice that the first condition 
also satisfies eq. (|19|) . and then it will result in a compact spectrum, in contrast, as we will see in the following, the 
second condition will yield a hierarchical spectrum. Before dealing with these two cases in detail, let us remark that 
without loss of generality, it is convenient to work in the basis where the RH neutrino mass matrix Mr is diagonal. 
Since Mr is symmetric, it can be brought to diagonal form M R lag = diag(Mi, M 2 , M 3 ) with real and positive entries 
by means of a unitary matrix W: 

M diag = w \ MrW * _ (25) 

In this basis we redefine the Dirac mass matrix as follows 

m D = m D W*. (26) 

From now on we will always work in this basis. 

1. Case 1: A 23 = 0, A 33 = 

In this case, we solve eq. (f2~T|) and expand the eigenvalues up to first order in 2^m. we obtain the following spectrum 
for the RH neutrinos 2 

Mi = 1^22^2, 

M 2 = \A 13 \m D im D3 , 

M 3 = \A 13 \m D1 m D3 . (27) 

With the reasonable assumption that |^4i3| and \A 22 \ are not very hierarchical, we see that it is possible to have 
Mi ~ M2.3 and, depending on the values of rn^ 119 renormalized at the leptogenesis scale, both mass orderings 
Mi < M2.3 or M2.3 < Mi are possible. In ref. pj only the Mi < M2.3 ordering was considered, a vanishing 6*13 was 
assumed, and both lepton flavour and heavier RH neutrino effects in leptogenesis had been ignored. Instead, as we 
will show, the ordering M2,3 < Mi can indeed occur, #13 7^ is a crucial condition to ensure the existence of solutions 
for the compact spectrum conditions, and as regards the lepton flavour and heavy RH neutrino effects, they must be 
included in order to obtain successful leptogenesis, and to guarantee that the result is reliable. 



2 Notice that the physical RH neutrino masses, eqs. I |27l . correspond to the absolute value of eigenvalues obtained from solving cq. I|21j l, 
Alternatively, one can also find the unitary matrix W in eq. ( 125ft which diagonalizes Mr up to first order in mpi . 



2. Case 2: A 2S = 0, A 22 = 

In this case, by proceeding as before up to first order in 2Lai ; we obtain the spectrum 

Mi = |A3 3 |m^ 3 , 

Mi = \A12\mDimD2, 

M 3 = \A 12 \m D1 m D2 . (28) 

Assuming also in this case that I-A12I and I.A33I are not exceedingly hierarchical implies M\ ^> M 2> s. Of course in this 
case, since the large contributions from mo 3 are not suppressed, we do not expect to obtain a compact RH spectrum. 
Nevertheless, in principle leptogenesis could still proceed at a scale M 23 <C Ajj thanks to the asymmetries generated 
in the decays of the two quasi degenerate states N 2 ^- Eventually however, we will find that in Case 2 on the one hand 
leptogenesis is unable to produce a sufficient baryon asymmetry, and on the other hand for the heaviest RH neutrino 
we always obtain M3 > 10 14 GeV which, under the requirement of perturbative Yukawa couplings, is in conflict with 
50(10) gauge coupling unification which instead suggests an intermediate vevs scale of order 10 11 GeV [34j (see 
however [35| for viable scenarios with an intermediate scale as high as 10 14 GeV). 

IV. RELATION WITH LOW ENERGY OBSERVABLES 

We have seen that by assuming conditions cq. (|20|) we have forcibly ended up with a quasi degenerate pair of RH 
eigenvalues. Before proceeding, let us stress that while setting the values of A 23 and A 23 to an exact zero has the 
virtue of simplifying the analysis, a generic compact RH neutrino spectrum can be obtained by fixing instead their 
values to any sufficiently small number as dictated by eq. (|18p. Doing this would lift the quasi degeneracy, but would 
still yield similar results. We will return to this point in Section |VJ One important point is that requiring that the 
matrix A satisfies some specific conditions gets reflected in specific relations between the low energy observables, and 
yields an enhanced level of predictability for the SO(10) model. Before studying which type of relations arise, it is 
useful to carry out a quick counting of the fundamental free parameters of the theory, and list the phenomcnological 
constraints that they should satisfy. The structure of the two symmetric matrices mp, and Mp is determined by 
two corresponding sets of fundamental Yukawa couplings between the fermion fields in the 16 and the Higgs fields 
respectively with vevs ~ Aew and ~ A#. This amounts to 6 + 6 + 1 = 13 real parameters corresponding to the two 
symmetric Yukawa matrices plus the ratio Apw/Ap which determines the seesaw suppression of neutrino masses, 
or in other words their absolute scale. Under the assumption of quark- lepton duality (see eq. [9]), the values of 
the 13 real parameters are constrained by the following observables in the up-type quark and neutrino sectors: the 
three quark masses m Ul m c , m f , the two neutrino mass-squared differences Am^, Am^, the three CKM mixing 
angles 8' 12 , 8' 23 , 8' 13 and the three PMNS mixing angles #12, #23, #13, which add up to a total of 11 constraints. Now, 
imposing on the complex elements of the matrix A two additional conditions, e.g. A 23 = A 33 = (or any other pair 
of conditions), implies that the set of 13 real fundamental parameters must satisfy two additional requirements, that 
in our case read Re(^23) = Re(A33) = 0. Thus the parameter space of the model remains completely determined 
allowing to obtain a quantitative prediction for the absolute neutrino mass scale m Vl . As regards the constraints 
on imaginary quantities, there are many fundamental complex phases, and only one measured observable, the CKM 
phase 8' . Nevertheless, as we will see, the structure of the conditions implies nontrivial relations between 8' and the 
three PMNS phases a, (3 and 8. 

In the following we assume a hierarchical and normally ordered spectrum for the light neutrino masses m*° 9 = 
diag(mi, m 2l m 3 ) with 

mi < m 2 < m 3 , (29) 

which is justified by the assumption of quark-lepton symmetry eq. ([§]). In the basis where charged lepton mass matrix 
is diagonal, the PMNS mixing matrix Upmns diagonalizes the effective neutrino mass matrix 

m v = U PMNS m u U PMNS . (30) 

We adopt for Upmns the standard parametrization in terms of 3 angles and three complex phases: 

C12C13 S12C13 si 3 e~ l5 

Upmns = [ -S12C23 - ci 2 s 23 si 3 e lS ci 2 c 23 - si 2 s 23 si 3 e lS s 23 a 3 ) x diag (l, e m , e l/5 ) . (31) 

S12S23 - ci 2 c 23 si 3 e tS -C12S23 - s 12 c 23 s 13 e lS c 23 c 13 



Here Cy = cos 9ij and Sij — sin 9ij, with i and j labeling families that are coupled through that angle (i,j = 1,2,3). 
Note that since the computation of the leptogenesis CP asymmetries involves several interfering amplitudes, the angles 
6ij cannot be restricted to the first quadrant, except for $13 that can be taken to be positive once the CP phase S is 
allowed to range between —it and w. 

According to the quark- lepton symmetry ansatz eq. ©, in the numerical analysis we take Vl = Vqkmi and 
accordingly we parametrize Vl with three angles and one phase, with a structure analogous to the first matrix on 
the right-hand-side (RHS) of eq. pip , distinguishing the angles and phase with a prime superscript: s' 12 , s 23 , s' 13 , S' 3 . 
The following discussion, however, is based on analytical expressions that get largely simplified by writing Vl in the 
approximate Cabibbo-like form: 



cos 9c sin 9c 
; ) - I 

' L 



VP = I -sm9 c cos 9 C I . (32) 



1 

We will then write our formulae in this approximation, keeping in mind however, that the full expressions have been 
used to obtain the numerical results. 

The matrix A eq. (fT6|) can be expressed in terms of the observables Vl , Upmns and mf, laa as 

A = (V L U PMNS T -^— (V L U PMNS y . (33) 



m 

AC) 



In the approximation Vl = V£ , the conditions A 23 = A 33 = yield the following two relations: 



m 2 _ 2ia _ (C12S23 + e. lS Si 2 s 1? ,C23) [s c c 12 si Z s 2 3 + e lS (s c s 12 c 23 - c e Ci 2 Ci 3 )] 

mi (-S12S23 + e J,5 Cl2Sl3C23) [-S c Sl2Sl 3 S23 + e" 5 (c c Sl 2 Cl 3 + S c Cl 2 C23)] 



!^3 e -20 _ 



C13C23 [s c Sl 2 C23 + C12 (-C c Cl3 + S c e lS Si 3 S 23 )] 



mi (-S12S23 + e i5 ci 2 si 3 C23) (s c ci 3 s 2 3 + c c e lS s 13 ) 

where s Cl c c = sm9c,cos9c- 

We see that by taking the absolute values of eqs. (l34l) - ([35|) we obtain two conditions that do not depend on the 
Majorana phases a, (3 (and are also even functions of 8 that depend only on its cosine). We can further eliminate 
m 2 and m 3 by using their relations with the solar and atmospheric mass-squared differences 4 : m 2 = m\ + Am, and 
m§ = m\ + Am^, obtaining: 

Am 

1 + — r = fi2([Oc,Oi2,023,Oi 3 ];cos6), (36) 



2 



/23([0C,012,023,013];COS(J). (37) 



The absolute neutrino mass scale mi appearing on the LHS of these equations represents the first unknown. On the 
RHS, /12 and /23 are two known (although non-transparent) functions of known mixing angles (that are listed within 
the squared brackets) and of the cosine of the second unknown, that is the Dirac phase S. These two equations might 
or might not have physical solutions (for example, given that Am^ > 0, in case the RHS of the first equation remains 
< 1 for all values of 5, there are no physically acceptable solutions). If solutions exists, these will corresponds to 
specific values of mi and of i5 with uncertainties determined by the experimental errors on the mixing angles. 
As regards the conditions on the complex arguments of eqs. (1M1) - (J3"5)) . they have the form 

a = gi 2 ([9 C : #12, 23 , #13, Awfl; 6, mi), (38) 

= gm([8c,8i2,Om,9i3,Aml]l6,m 1 ), (39) 



3 Taking Vl ~ VcKM implies that the large leptonic mixing observed in the low-energy sector should be a consequence of a see-saw 
enhancement of lepton mixing. Such an enhancement requires a strong (quadratic) mass hierarchy of RH neutrinos, or a Mr structure 
with large off-diagonal entries 13611 . 

4 The relation with the atmospheric mass-squared difference rrig — c 2 s rn^ — sj.m\ = Am^ (see ref. Q) can be equally well used with 
irrelevant numerical differences. 



Quark sector 


Neutrino sector 


m„(A) 


0.00067 GeV 


Amir (A) 


11.86 x 10" 5 eV 2 


TOe(A) 


0.327 GeV 


Ami! (A) 


3.84 x 10" 3 eV 2 


m t (A) 


99.1 GeV 






®\2 


13.02° 


#12 


34.4° 


823 


2.35° 


823 


42.8° 


#13 


0.20° 


013 


5.6° 


8' 


1.20 rad 







TABLE I: Input parameters. We use the up-quark masses renormalized to the scale A = 10 9 given in Table IV in Ref. [371 ]. 
Neutrinos mass squared differences are taken from the global fit in Ref. [15[ and renormalized to the scale A with a multiplicative 
factor r 2 with r — 1.25 according to the prescription in Ref. 13811 . The CKM mixing angles 6'^ and CKM phase 8' are derived 
from the values of the Wolfenstein parameters given in Ref. [23|. The PMNS mixing angles are taken from the global fit in 
Ref. [13|. Renormalization effects for the CKM and PMNS parameters have been neglected. 
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FIG. 1: Plots of mi as a function of 8 according to eq. (|36|) (solid blue line) and eq. (|37p (dashed violet line). The points 
of intersection represent pairs of possible solutions for (m vi ,8). Left panel (+|#i2|, — 1 6*23 1 ) , right panel (— |0ia|, — |023|), with 
(|6>i 2 |,|6>23|) = (34.4°, 42.8°). 



where c/12 and g±3 are again known functions, and thus a and ft can be determined in terms of mi, 6 and of the known 
mixing angles and mass squared differences. This completes the determination of all the low energy observables in 
terms of measured quantities, and through eq. (|15[) also fixes the RH neutrinos mass spectrum. 

To be more precise, given that the signs of #12, #23 and #13 are not determined in oscillation experiments, depending 
on the possible choices ±6^ the two eqs. (|36p - ([37)l represent in principle 2 3 = 8 conditions. However, the PMNS phase 



always appears together with 6*13 in the combination Si3e = —Sis^' so that if S is a solution for +#13, S ± tt is 
a physically indistinguishable solution for —#13, and this reduces the eight possible pairs of equations to just four. 

Once the simplification Vl —$■ V£ is dropped, due to the presence of the CKM phase 5' eqs. (j36]) - (l37T) acquire 
a (mild) dependence also on sin 5, meaning that for each one of the four possibilities (±#12, ±#23) we can have two 
nonequivalent solutions corresponding to values of S of opposite signs. An example of this situation is illustrated in 
figEQfor the two cases (+,-) = (+|6>i 2 |, -|6» 23 |) and (-,-) = (-|6»i 2 |, -|6» 23 |). The two relations eq. K35J} and eq. ([57)1 
correspond to two different curves m\{8) that are plotted respectively with the solid blue lines and the dashed violet 
lines and intersect in two points that are the solutions to the system of constraints. Notice that a solution to these 
constraints does not always exists i.e. when the two curves do not intersect. This happens for example in some of the 
scenarios in Case 2, as shown in Table ITiTl that has therefore less entries than Table [TT1 of Case 1. 

The input parameters of our numerical analysis are listed in Table U] For the eigenvalues of mo we use the values 
of the up-quark masses renormalized to the scale A = 10 9 GeV (~ Mr), given in Table IV in Ref. [37]. The relevant 
values of Mr we find are given in Table |TTJ Neutrinos mass square differences are taken from the global fit in Ref. 11511 
and renormalized to the scale A with a multiplicative factor r 2 with r = 1.25 according to the prescription in Ref. [38l |. 
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FIG. 2: Left: the iVi CP asymmetries t\ e (solid blue lines) ei M (dashed violet lines) ei T (dotted green lines) for the solution 
labeled ( — , — ) in Table lU (Case 1) for positive (thick lines) and negative (thin lines) values of the Dirac phase S, as a function 
of the Vl phase S' . Right: the Ni electron flavour washout projector P\ e for positive (thick line) and negative (thin line) values 
of 5 as a function of 5'. The vertical dot-dashed lines correspond to the value <5' = 5c km — 1.20 rad used in the numerical 
analysis. 



The CKM mixing angles 6L and CKM phase 5' are derived from the values of the Wolfenstein parameters given in 
Ref. [39J , renormalization effects for these angles are small and have been neglected. 

As regards the PMNS mixing angles, recent fits to oscillation neutrino data suggest a small but nonvanishing value 
for #13. In our scenario, having #13 7^ is of fundamental importance because only under this condition the Dirac 
phase 8 will enter the constraining equations eq. (l34l) and eq. (|35|) . providing enough free parameters to allow for 
a numerical solution, so let us discuss this specific quantity a bit more in detail. With the assumption of normal 
ordering mi < m 2 < m 3 , and with 3tr errors for the three-flavour neutrino oscillation parameters, the following results 
have been reported: 



sin 2 6 13 = 0.013+^^[ll, sin 2 13 = 0.025±g;g| 5 [I3], sin 2 13 < 0.043 [H 



(40) 



The last result is a 3cr upper limit estimated in the framework of the so called GS98 solar model with the Ga capture 
cross-section of Ref. 0. At la the same data give sin 2 $ 13 = 0.0095to:o(r7 corresponding to 6 13 = (5.6t|^)° [ll]. We 
use such best fit value, and also for the other two angles we adopt the results of the global fit in Ref. fla ]. 

Our results for the possible values of the still unmeasured low energy parameters mi, 8, a and j3, and for the RH 
neutrino masses evaluated according to eq. (fT5|) , are collected in Table |TT] and in Table Mil Here onwards we always 
arrange the ordering of RH neutrino masses according to Mi < M 2 < M 3 . The first four lines of Table HT1 list the 4+4 
possible solutions for the conditions of Case 1 A 23 = A 33 = 0. In each line we list the two solutions corresponding 
to positive and negative values of 8. Note that the numerical differences between the absolute values of each pair of 
solutions for 8 are small, since they correspond to effects suppressed by s' 13 . However, even such a small difference 
can have a non negligible impact on the value of the leptogenesis CP asymmetries. Real parameters like mi and Mj 
also come in pairs with very close values, but in this case the differences are numerically irrelevant so that a single 
approximate value is displayed. This situation is illustrated in Fig. [2j the left panel depicts the three N\ flavoured 
CP asymmetries (see the next section) €\ e (solid blue lines) ei p (dashed violet lines) e\ T (dotted green lines) for the 
solution labeled (— ,— ) in Table [TTI (Case 1) for positive (thick lines) and negative (thin lines) values of the Dirac 
phase S, as a function of the Vl phase 8' ■ We see that for the CP asymmetries which are very sensitive to the values 
of the complex phases, the two different solutions for 8 induce very large effects, for example they swap completely 
the signs of e\^ and e\ T . In the right panel, as an example of one important real leptogenesis parameter, we have 
plotted the electron flavour washout projector P\ e for the RH neutrino N\ (see next section), for positive (thick line) 
and negative (thin line) values of 8 as a function of 8' . We see that in this case numerical differences are irrelevant. 
In both panels the vertical lines correspond to the value 8' = Sckm = 1.20 rad that we have used in the numerical 
analysis of Tables HHHED 

The last line in Table [TTI labeled with (—,—)* gives the results obtained for that case when A 23 , A 33 are set to small 
but nonvanishing values, that we have (arbitrarily) chosen as (cfr. eq. (|18[) ): 
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(#12, #23) 


5 


mi (10^ 3 eV) 


a 


P 


(Mi, M 2 , M 3 ) (10 9 


GeV) 


(+,+) 


(1.43,-1.46) 


9.7 


(-1.47,1.46) 


(-0.18,0.19) 


(3.5, 3.9, 3.9) 


(-,+) 


(2.60, -2.63) 


2.9 


(1.52,-1.51) 


(1.25,-1.22) 


(3.0, 8.7, 8.7) 


(+,") 


(1.91,-1.89) 


8.2 


(1.46,-1.44) 


(0.13,-0.14) 


(4.0, 4.0, 4.3) 


(-,-) 


(0.74,-0.71) 


3.1 


(-1.49,1.48) 


(-1.19,1.17) 


(3.5, 7.9, 7.9) 


(",-)• 


(0.74,-0.71) 


3.1 


(-1.49,1.48) 


(-1.19,1.17) 


(3.5, 7.8, 8.0) 



TABLE II: The eight possible solutions to the RH neutrino compact spectrum conditions of Case 1: A23 = A33 = 0. The real 
parameters mi and Mi also come in pairs as the complex phases 5, a and j3. A single approximate value is displayed because 
the difference between the two values is numerically irrelevant. For the (+, — ) solution, the almost degenerate RH neutrinos 
are the lighter ones, while in all the other cases are the heavier ones. In the last line, labeled with an asterisk ( — , — )*, we give 
the results for the (— , — ) solution in which A22, and ^33 are set to arbitrary small, but nonvanishing values (see text). 



We see that the changes in S, mi, a and f3 with respect to the A23 — A33 — case in the fourth line remain below 
the precision of the table, and in any case are way too small to be seen experimentally. The same happens for all the 
other cases, and we can thus conclude that the predictions obtained with our simplified conditions eq. ([20)1 hold for 
each class of compact spectrum solutions. The last column of this line however, makes apparent that the resulting RH 
spectrum is just compact rather than degenerate, which justifies talking about a third way to SO(10) leptogenesis. 

The conditions A23 = A22 = of Case 2 have only 2+2 solutions, that are listed in Table IIIII In this case there 
is a very large hierarchy O(10 7 ) between the two almost degenerate RH neutrino masses Mi, 2 and M3 so that, as 
expected, the RH spectrum is not compact. Perturbativity of the Yukawa couplings then implies Kr <^ 10 14 GeV. 
Accommodating such a large intermediate scale might be problematic in the SO(10) model. Furthermore, as we will 
see in the next section, in this case leptogenesis is unsuccessful for both the (— , +) and (— , — ) solutions. 

To conclude this section, we have seen that by forcing the SO(10) model to produce a compact RH neutrino 
spectrum we obtain a scenario in which all the parameters relevant for leptogenesis remain determined in terms of the 
set of low energy observables that have been already measured. Most remarkably, as we will see, besides predicting 
values for the absolute neutrino mass scale mi, the CP violating phases, and the RH neutrino spectrum, the size and 
the signs of the flavoured CP asymmetries relevant for leptogenesis are also fixed, and allow to predict the size and 
sign of the cosmological baryon asymmetry generated through leptogenesis. Such a level of predictability is indeed 
quite unusual in see-saw inspired scenarios. Clearly, verifying if the baryon asymmetry yield of leptogenesis is in 
agreement with observations will now represent the major test of our scenario. This is the task that we are going to 
address in the next section. 



(#12, #23) 


s 


mi(10- 3 eV) 


a 


P 


Mi, 2 (10 7 GeV) 


M 3 (10 14 GeV) 


(-+) 


(2.66, -2.69) 


2.0 


(-1.51,1.52) 


(-0.21,0.23) 


3.1 


3.6 


(-,-) 


(0.35, -0.32) 


2.0 


(1.53,-1.54) 


(0.19,-0.21) 


3.3 


3.0 



TABLE III: Same than Table (TT) for the four possible solutions to the conditions of Case 2: A23 = A22 — 0. 



LEPTOGENESIS 



In the basis for the Dirac mass matrix as in eq. (l26|) . the CP asymmetry in the decay of the RH neutrino N{ 
(i = 1, 2, 3) to a lepton £ a (a = e, fj,, t) is given by [41[ 
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where v — 174 GeV is the EW vev and 5 

9{x) 



1-x 



(i-*) 2 + (w-^)' 

1 -x 



+ 1- (l + x)lo{ 



1 + x 



a-*) a + (ft-*£t)' 



(43) 



are loop functions with I\ 






{m D rho)u the total iVj width. The expression in the second line of eq. (|4"2"j) corre- 
sponds to the lepton-fiavour-violating but lepton-number-conserving self-energy diagram. It vanishes when summed 
over a so it does not contribute in the one-flavour approximation, but plays an important role [441 ] when, as in our 
case, leptogenesis occurs in the flavoured regime. In eqs. (|4"3"|). g(x) and the first term in the square bracket of f(x) 
come from the self-energy contributions with the resonant condition given by 



± 



M, 



'Aft 



(44) 



while the remaining terms in f(x) correspond to contributions from the vertex diagram. The resonant condition 
cq. (|4~4")) gives 



f(x) ~ y/xg(x) 



\Mi ^ M k 



(45) 



where in f(x) we have ignored the subleading contributions of the vertex diagram. However, in our case although the 
degeneracy conditions M^ ~ M% or M± ~ M.% are approximately fulfilled, we never reach a fully resonant regime as 

defined by eq. (|4"4"jl . and hence ignoring the "regulator" term I jf- — Xj-jf- ) in eqs. (|43| only yields negligible numerical 

differences. 

In order to calculate the baryon asymmetry, we need to solve a set of Boltzmann equations (BE) (we refer to [6j 
and references therein for details). By including for simplicity only decays and inverse decays, the BE for the RH 
neutrino densities Y^ i and for Y Aa , that is the asymmetry density of the charge B/3 — L a normalized to the entropy 
density s, can be written as: 
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1 N 
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2 V Y e eq 



veq 
1 H 



(46) 



where Y^ ? 



45 



z 1 K,2{z) is the equilibrium density for the RH neutrinos with g* = 106.75 and K-2 the second order 
modified Bessel function of the second kind, 2Y^ q = Y^ q = 47r 1 2 5 are respectively the equilibrium densities for lepton 
doublets and for the Higgs, and the integration variable is z — M/T with T the temperature of the thermal bath, and 
M = M\ t 2,3 the mass of the decaying neutrino. In the above, we have defined Ya q = ^as/3 — ^al q with Yal q the 
total lepton density asymmetry in the a flavour which also includes the asymmetries in the RH lepton singlets. Since 
RH neutrinos only interact with lepton doublets, the RHS of the second equation of eqs. (J46]) involves only the LH 
lepton doublets density asymmetry in a given flavour a, Y A e a = A a pY Aa with A a p the flavour mixing matrix 8] given 
in eq. (Unj) below. In the same equation, we also define Yah = CpY Afj the Higgs density asymmetry with Cp [45[ also 
given in eq. (|4T))) and jN ia = PialNi (no sum over i) where Pi a projects the decay rate over the a flavour, that is, it 
corresponds to the branching ratio for Ni decaying to £ a , and can be written as 



P; r 



_\ / la. 

v / a 



(47) 



For the resonant terms, we used the expressions from Refs. 1421. 143 
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Let us also introduce the rescaled decay width 

M? Mi ' v ' 

which is also known as the effective washout parameter, that parametrizes conveniently the departure from thermal 
equilibrium of A^-related processes (the larger m„ the closer to thermal equilibrium the decays and inverse decays 
of Ni occur, thus suppressing the final lepton asymmetry). Finally, the combination Pi a fhi projects the washout 
parameter over a particular flavour direction, and determines how strongly the lepton asymmetry of flavour a is 
washed out. 

Leptogenesis becomes possible when the thermal bath temperature approaches the value of the mass of the decaying 
RH neutrino, that becomes non relativistic and can decay. However, when the washout parameter is large m <^ 10 , 
at T ~ Mi the RH neutrinos are actually in equilibrium and no asymmetry can be generated. In our case the RH 
neutrinos are coupled rather strongly to the thermal bath (m 3> 10~ 3 see e.g. Table HV|) and in this case the generation 
of the bulk of the lepton asymmetry is delayed down to much lower temperatures: for fh ~ 6 x 10 -2 for example one 
can estimate z = M/T ~ 8 |46j). Thus, the range of temperatures where the lepton asymmetry is generated falls 
well below 10 9 GeV, where both the r and /x Yukawa interactions are presumably in equilibrium jg|. In this regime 
all the three lepton flavours are then distinguished, and their dynamical evolution must be followed separately. The 
A flavour mixing matrix and the C vectors allow to accomplish this task, and in our temperature regime are given 
by0 

-, / -906 120 120 

A = 75 -688 28 



2148 



75 28 -688 



C = --^(37,52,52). (49) 

Once the final asymmetries in the lepton flavour charge densities Ya„ are obtained by solving numerically the 
BE eq. (|4"6")l . the baryon asymmetry generated through leptogenesis is given by [47j 

^AB = ^E F A Q - (50) 

a 

The resulting prediction should then be confronted with the experimental number. The most precise experimental 
determination of Yab is obtained from measurements of the cosmic microwave background (CMB) anisotropics. A 
fit to the most recent observations (WMAP7 data only, assuming a ACDM model with a scale-free power spectrum 
for the primordial density fluctuations) [48j], when translated in terms of Yab gives at 95% c.l. [49j 

Yab B = ( 8 - 79 ± °- 44 ) x 1(r11 - ( 51 ) 

A. Numerical Results 

For all the phenomenologically viable cases stemming out from our scenario, the complete set of high energy 
parameters required for computing the baryon asymmetry yield of leptogenesis is predicted. The RH neutrino masses 
are listed for Case 1 in the last column in Tabic [Hi and for Case 2 in the last column of Table IIIII (of course, the 
precise numerical values of the masses, rather then the approximate values listed in the tables, are used for the 
leptogenesis computation). As regards the flavoured CP asymmetries, they are computed according to eq. (|42)l and 
the corresponding results for our two cases are listed in Table IIVI and [Vj Case 1 (Table IIVI) results in a very compact 
spectrum of RH neutrinos (and with a pair of almost degenerate states if the exact conditions A23 — A33 = are 
imposed). For solutions (+, +) and (+, — ) the largest mass differences remain at the 10% level, while for solutions 
(— ,+) and (— ,— ) they reach a factor of a few. Under this conditions it is mandatory to include the contributions 
from the heavier RH neutrinos, that in the first case (of tiny mass differences) can affect the results for a factor up to 
O(10 3 ). For Case 1 ( Table HV]) . that contains sub-cases in which leptogenesis can be successful, we give the complete 
list of the Ni flavoured parameters: the CP asymmetries for the three flavours are given in columns 2-4, the washout 
flavour projectors in columns 5-7, and the effective washout parameters in column 8. The final results obtained by 
integrating the BE and by converting Y& B _ L — ^2 a ^A Q into Ya b through eq. (|50|) are given in the last column of 
the two tables. Columns 2-4 in Table IIVI show that asymmetries of different signs are produced in different lepton 
flavours (an example of this situation is also depicted in the left panel in Fig. [5]), while columns 5-7 show that a certain 



(#12, #23) 


ei e (lin 5 ) 


cim(1CT 5 ) 


ei T (10~ 5 ) 


Pie 


P* 


Pit 


mi(10" 2 eV) 


Y AB (W- 10 ) 


(+,+) 


(0.9,0.5) 


(3.4,-5.0) 


(-4.3,4.5) 


10~ 3 


0.03 


0.97 


74 


(-0.3,-0.08) 


(-+) 


(0.041,0.083) 


(-0.52,0.41) 


(0.47, -0.49) 


0.03 


0.50 


0.47 


6.2 


(-0.8,-1.5) 


(+,") 


(-0.001,-0.001) 


(0.003, -0.001) 


(-2.5,2.6) 


10" 4 


io- 3 


1.0 


IO 5 


(0.03,0.07) 


(--) 


(-0.16,-0.06) 


(-0.89,1.14) 


(1.05,-1.09) 


0.03 


0.47 


0.50 


6.7 


(2.3,0.32) 


(",-)* 


(-0.05,-0.04) 


(-0.05,0.15) 


(-0.15,0.15) 


10~ 3 


0.06 


0.94 


70 


(0.43,0.38) 



TABLE IV: The leptogenesis parameters corresponding to the different solutions of Case 1. 
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hierarchy exists between the flavoured washout parameters. Given the highly non uniform pattern in flavour space of 
the relevant leptogenesis quantities it is clear that no analytical expression based on the single flavour approximation 
would produce a reliable result, and we can firmly conclude that lepton flavour dynamics is of crucial importance for 
studying leptogenesis in the SO(10) model 6 . 

For Case 2, given that all the solutions consistent with the low energy constraints eventually fail the leptogenesis 
test, we just give in Table IV1 the total CP asymmetries and the approximate values of the total washout parameters 
for the two quasi degenerate RH neutrinos iVi.2- This reduced set of figures is however sufficient to conclude at a first 
glance that with CP asymmetries of 0(1CP 8 ) and washout parameters of 0(1 eV), no flavour dynamics could rescue 
leptogenesis from a quantitative failure. 

Our results are collected in the last columns of Tables IIVI and Table [V] Althou gh w e are using somewhat simplified 
BE in which thermal corrections [38] , scatterings and CP violation in scatterings (50H52J ] , and other subleading effects 
are neglected, the estimates of the final baryon asymmetry we obtain should be sufficiently accurate for our scopes. 
For example, we have checked that including scatterings and CP violations in scatterings introduces a < 25 % effect, 
which is by no means crucial to test the scenario. There are in fact other important sources of uncertainties: in our 
analysis we are using fixed central vales for all the input parameters, and it goes without saying that the final value of 
Yab will be affected by the experimental uncertainties. Even more importantly, there are also theoretical uncertainties 
stemming from deviations from the exact quark- lepton symmetry ansatz eq. Q, as well as from deviations from the 
exact zeroes in the conditions A 2 3 = A 33 = 0, which are obviously difficult to quantify. Therefore, we will be contented 
to require that a successful prediction of the BAU, besides having the correct sign, should approach the experimental 
result eq. ([5"Tj) only within a factor of a few. For Case 1, we obtain four solutions with the wrong (negative) sign of the 
BAU, and other four with the correct sign. They are listed in Table ITVl However, only the two solutions in the fourth 
line of the Table are sufficiently close to the experimental value eq. (f5"Tj) to be all phenomenologically acceptable. For 
this two solutions we give in the last line of the Table the values of the leptogenesis parameters for Ni , and of the final 
baryon asymmetry when the exact zeroes in eq. (|20|) are lifted to small but not vanishing values as given in eq. (|4T[) . 
We see that although the final value of Yab is sensitive to this change, it still remains within a factor of two from the 
measured central value eq. (|5T|) . 

As we have already said, in this SO(10) scenario the leptogenesis efficiency gets largely enhanced by flavour effects, 
and it is then worth asking what would happen if the bulk of the lepton asymmetry, rather than in the three flavour 
regime, is generated when only the r Yukawa coupling mediated in-equilibrium reactions, and the number of relevant 
flavour is reduced to two. In the two flavour and strong washout case, estimating Yab is more subtle because there can 
be protected directions in which the asymmetry generated by ^2,3 is not erased by N\ washouts [HIlLII]- We follow 
a simplified approach that neglects this phenomena, and thus gives a conservative estimate of the final asymmetry, 
obtaining for the (— ,— ) case Yab ~ (0.3,-0.3) x 10~ 10 . In the first case the asymmetry gets reduced, but still 
remains within a factor of three from the experimental number; however, in the second case the asymmetry changes 
sign, which again shows the importance of a proper treatment of flavour dynamics. 

As regards Case 2, we see that the two right-sign solutions yield a baryon asymmetry that is too small by almost 
three order of magnitudes. This suppression of Yab is due to two different reasons: firstly the CP asymmetries are 
exceedingly small because the imaginary parts of the relevant combinations of couplings are strongly suppressed, and 
secondly the washout parameters are rather large, and imply an almost in-equilibrium dynamics that impedes building 
up any sizable density asymmetries. 



For example, an analysis in some aspects similar to ours, but in which flavour dynamics is neglected, has been carried out in Ref. 1221 . 
Their 'Special case III' is similar to our Case 1, but the conclusions are opposite. This is most likely due to the enhancements of the 
leptogenesis efficiency from flavour effects. 
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(#12, ^23) 


ei (1(T 8 ) 


e 2 (W~ S ) 


tii, 2 (eV) 


Y AB (10- 14 ) 


(-,+) 


(5.4,-5.5) 


(5.4,-5.5) 


1.7 


(-1.8,1-8) 


(-,") 


(-4.3,5.0) 


(-4.3,5.0) 


1.6 


(1.5,-1.7) 



TABLE V: The leptogenesis parameters corresponding to the different solutions of Case 2. 



In conclusion, the SO(10) model constrained by the assumption of the quark-lepton symmetry in eq. © and by 
the compact RH neutrino spectrum conditions in eq. (|20[) . when confronted with the results from neutrino oscillation 
experiments, and with the requirement of successful leptogenesis, yields predictions for the yet unknown low energy 
neutrino parameters, that are summarized in the following two possibilities 



(— ) 



3x lCT J eV, 



±0.7, 



Tl.5 ; 



/?- Tl-2, 



(52) 



which correspond to either the upper or lower sign of the three phases. 

With the numerical results listed in eq. (f52"j) another low energy observable can be predicted, that is the neutrinolcss 
double beta decay effective parameter 



m,, 



(m„)n = ^{Up MNS )u{m v )f aa {ul 



PMNSl* 1 ' 



(53) 



for which we obtain \m ee \ 5i 2 x 10 3 eV that, as could have been expected for hierarchical and normal ordered 
neutrino masses, remains well below the sensitivity of all ongoing and planned experiments |53j |. 



VI. DISCUSSION AND CONCLUSIONS 



The predictive power of our scenario is spelled out in clear in eq. (j52|) . and such a high level of predictability calls 
for an explanation. The crucial point is that in our study there are no free parameters: everything is fixed in terms 
of the low energy neutrino observables and by the additional assumption of quark-lepton symmetry cq. §§§ and by 
the compact RH spectrum conditions A23 — A33 = 0. The only freedom left over by these latter constraints is a 
discrete one, and corresponds to the signs of the two angles #12 and #23, for each choice of which there are in turn two 
solutions, corresponding to positive and negative values of the phase S. Given that there is no free parameter that can 
be adapted to fit the observed value of the BAU, we find intriguing that among the discrete set of eight possibilities 
of Case 1, in two cases the leptogenesis yield of baryon asymmetry is in acceptable agreement with observations. 

To summarize the main results of the paper, we have first shown that in the SO(10) seesaw model it is technically 
possible to arrange for a compact RH neutrino spectrum, and this in spite of the fact that the SO(10) neutrino Dirac 
mass matrix is characterized by a hierarchy between its eigenvalues that is much stronger than the one observed for 
the the light neutrinos, a situation that would naturally call for a compensating large hierarchy in the RH masses. 
We have argued that this possibility can be implemented in a consistent way only if the PMNS mixing angle #13 is 
nonvanishing, since only in this case we have at disposal the Dirac phase 5 as an additional free physical parameter 
that can cope with satisfying the compact RH spectrum conditions. The counting of free parameters is a subtle point: 
clearly our construction relies quantitatively on the assumption of a strict quark-lepton symmetry mjj — m u and 
Vl = VcKMt and one might argue that even in the absence of 5 one could be able to find solutions by modifying 
these assumptions. Nevertheless, in SO(10) ttid and Vl are intrinsically related to m u and Vckm, simply by the fact 
that all fermions of each family, including the RH neutrino, are assigned to the same irreducible representation of 
the group. The important point is that while the precise form of the quark-lepton duality relations can be changed 
according to the amount of (family dependent) contamination in the EW breaking sector from 10 and 126 vevs, 
still mu and Vl cannot be regarded as independent from the corresponding quantities in quark sector, since for any 
fixed pattern of vevs, some specific relation between them remains fixed. Therefore, the compact spectrum conditions 
together with any specific assumption about lepton-quark Yukawa relations, always yields a scenario where the values 
of the yet unknown neutrino parameters can be predicted directly in terms of known quantities. This is not so for 
the absolute neutrino mass scale mi which is a true free parameter, since it is essentially determined by the ratio 
Aew/A-r where the scale A# is free, nor for the PMNS phase 5 since the complex phases in Mr, that are unrelated 
to Vckm and thus are also free, concur to determine its value. 

As regards the compact RH spectrum conditions A23 = A33 — eq. (|20l) . they should be understood with a grain 
of salt. Rather than corresponding to exact zeroes, the assumption is that to a good approximation the values of 
these entries are negligible. This is spelled out in eq. (fTHj) and eq. (flT))) . We have also tested the effects of lifting 
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the exact zeros to the small nonvanishing values given in eq. (|4"Tj) finding that (i) the quasi degeneracy in the RH 
neutrino mass eigenvalues is removed, resulting in a generic compact spectrum; (ii) the predictions for the measurable 
low energy parameters (the absolute neutrino mass scale and the CP violating phases) are not changed; (iii) the 
effects on the final value of Tab remain under control. We stress again that we have not put forth any theoretical 
explanation, as for example a symmetry argument, for why the entries A23 and A33 should be particularly suppressed, 
and this implies that the fact that the four (real plus imaginary) conditions can be fulfilled only if specific quantitative 
relations between mi, S, a and /3 are satisfied, should not be confused with a parametric (i.e. functional) dependence 
like S = 5(77*1) or a = a(mi,5), which is something that our SO(10) scenario certainly does not give, but should 
rather be regarded as numerical accidents. 

It is likely that in the not too far future the values of mi and of S will eventually be measured, and therefore 
the specific scenario we have been exploring, and whose predictions are summarized in eq. (|52l) . is straightforwardly 
falsifiablc. Of course, by modifying the form of the quark-lepton symmetry relations one would obtain numerically 
different predictions. However, any different assumption would result in the same level of predictability, and in 
particular it will have to pass the leptogenesis test which, as we have seen, is a highly nontrivial requirement. It is 
certainly conceivable a situation in which no assumption will be able to reproduce the measured values of mi and 8 
while simultaneously pass the leptogenesis test. We can then conclude that the leptogenesis scenario based on SO(10) 
with a compact RH neutrino spectrum is a testable physical hypothesis. 

Note added 

After this paper was published in the arXiv . org database [54[, the Daya Bay reactor neutrino experiment announced 
the measurement of a non-zero value for the neutrino mixing angle #13 with a significance of 5.2 standard deviations: 
sin 2 26»i3 = 0.092 ± 0.016(stat.) ± 0.005(syst.) [55J]. I n A P ri l 2012, the RENO experiment also reported a non zero 
value sin 2 26*13 = 0.113 ± 0.013(stat.) ± 0.019(syst.) [56| consistent with the Day Bay result. As we have explained 
in the paragraph above eq. (1401) . #13 ^= is a mandatory condition for the consistency of our scenario, which is now 
ensured by the Daya Bay and RENO results. The experimental central values #13 = 8.8° (Daya Bay) and #13 = 9.8° 
(RENO) are larger than our reference value #13 = 5.6° [15[ (see below eq. (|4T)|) ). and this will slightly change the 
numbers in eq. (|52|) . However, the conclusions of the leptogenesis analysis that are based on the sign of the baryon 
asymmetry and on its value only within a factor of a few, will not be changed. To give an example, for #13 = 8.8° we 
obtain for the first of the two cases labeled (-, -) in Table |TV]Tab = 2.9 x 10 -10 instead than Y AB = 2.3 x 10~ 10 
which clearly implies the same conclusions. 
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